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Universal logarithmic terms in the entanglement entropy appear at quantum critical points 
(QCPs) in one dimension (ID) and have been predicted in 2D at QCPs described by 2D con- 
formal field theories. The entanglement entropy in a strip geometry at such QCPs can be obtained 
via the "Shannon entropy" of a ID spin chain with open boundary conditions. The Shannon entropy 
of the XXZ chain is found to have a logarithmic term that implies, for the QCP of the square-lattice 
quantum dimer model, a logarithm with universal coefficient ±0.25. However, the logarithm in the 
Shannon entropy of the transverse-field Ising model, which corresponds to entanglement in the 2D 
Ising conformal QCP, is found to have a singular dependence on replica or Renyi index resulting 
from fiows to different boundary conditions at the entanglement cut. 
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The use of quantum information concepts to under- 
stand many-particle systems near a quantum critical 
point (QCP) has grown rapidly since the 1994 calcula- 
tion of entanglement entropy in one-dimensional (ID) 
critical systems pj. Entanglement entropy shows a uni- 
versal logarithmic divergence at ID critical points de- 
scribed by 2D conformal field theories (CFTs) and at 
infinite-randomness ID critical points Recent work 
has studied generalizations to Renyi entropy and the en- 
tanglement spectrum [1] . Terms proportional to log L for 
a subsystem of size L in an infinite background are es- 
pecially important as their coefficients are independent 
of the microscopic lattice spacing and hence potentially 
universal. Entanglement can be used to develop a classi- 
fication of ID interacting quantum systems and deter- 
mines the difficulty of numerical simulation of ID QCPs 
via matrix product state methods [5]. 

Above one spatial dimension, there are few general 
results on critical entanglement entropy. For critical 
points described by a c? -I- 1-dimensional CFT there is 
a conjecture for some geometries from the AdS/CFT 
correspondence fZl for the powers of length that appear 
in the entanglement entropy: logarithmic terms appear 
generically only in odd spatial dimensions. This conjec- 
ture agrees with ID results, as well as perturbation the- 
ory [5] and variational methods |1] above ID. For another 
class of systems, "conformal quantum critical points" 
(CQCPs) in d = 2 [TO], whose ground state wavefunc- 
tions are related to 2D rather than 3D CFT's, nonper- 
turbative analytical and numerical results are possible 
and are the focus of the present work. An example of 
a CQCP is the critical state of the square lattice quan- 
tum dimer model originally introduced as a model 
for high-temperature superconductors. 

We investigate logarithmic terms in the von Neumann 
entropy S and Renyi entropies 5„ for CQCPs associated 
with the compact boson (which has CFT central charge 
c = 1) and the Ising model (c = 1/2). One result is 



the first explicit observation of a universal entanglement 
entropy logarithm above ID. For the c = 1 boson, the 
results confirm a conjecture from the initial study of en- 
tanglement at CQCPs [H], which derived a formula for 
a universal, geometry-dependent logarithmic correction 
to the leading area law. Recent work [TBHT5] on order- 
unity (0(1)) terms called into question the validity of 
that formula as for the 0(1) terms there are subtle issues 
related to compactification that took several years and 
several groups to resolve. We give an improved derivation 
in supplementary material explaining why subtle factors 
affecting the 0(1) terms are irrelevant for the logarithm. 

However, for the Ising CFT we find that the coeffi- 
cient of the logarithm is discontinuous as the Renyi in- 
dex n passes through n = I. We explain this behav- 
ior analytically by arguing that the correct "defect line" 
used in the calculation of the entropy changes discon- 
tinuously with n and is not the combination of one free 
field and n — 1 Dirichlet fields as conjectured in Ref. [12] 
and confirmed for the boson. A similar but less singular 
discontinuity in n was previously found numerically to 
exist for the 0(1) term [161. Our numerical approach to 
logarithmic terms uses large-scale Time-Evolving Block 
Decimation (TEBD) [T7] calculations to implement the 
same mapping used there between entanglement entropy 
at 2D CQCPs and the Shannon entropy of ID spin chains. 
The compact boson case relevant to the quantum dimer 
model seems to be well understood, but the results here 
demonstrate that the correct defect boundary condition 
can be complicated and n-dependent. 

We first define some basic concepts. The bipartite 
entanglement spectrum {pa} of a pure quantum state 
consists of the eigenvalues of the reduced density matrix 
of either subsystem. The von Neumann entropy S and 
Renyi entropies Sn are defined by 

^^-^p.logp,, 5„ = ^log(^p^). (1) 



2 



(a) 2+lD CQCP 



A 



D 



(b) 2D CFT 



T - 

(c) ID spin chain 
Pa = \{(T\0)f 

CKMXD-O-O 



o-o-o-o-o-o 



FIG. 1. a) The entanglement spectrum of the CQCP, Po-, 
in a strip geometry is in one-to-one correspondence with field 
configurations {a} on the entanglement cut. These probabil- 
ities can be alternately understood as b), partition functions 
for a constrained CFT in the strip geometry, or c), ground 
state amplitudes of a quantum chain. 



The entanglement spectrum {pa} of a CQCP ground 
state in strip or cylinder geometries can be understood 
both in terms of the associated 2D CFT or as proba- 
bility amplitudes of a l-flD spin chain, as illustrated in 
Fig. [l] The entanglement spectrum {po-} is in one-to- 
one correspondence with CFT field configurations {a} 
along the entanglement cut. If the CQCP lies on a strip 
or cylinder, the 2D CFT naturally defines a ID quan- 
tum Hamiltonian H via the transfer matrix formalism. 
'Time' is chosen parallel to the strip (cylinder), with the 
entanglement cut at a fixed time t — 0. For universal 
properties, we can take the spatial direction to lie on a 
lattice. In the limit of an infinite strip the ground state 
probabilities are identical to the entanglement spectrum 
of the associated CQCP [13]. The Shannon entropy of 
the spin chain, defined as in Eq. ([T]) with p„ — |(ct|0)P, 
is a basis dependent measure of disorder in the ground 
state. The Shannon entropy of the spin chain is then 
equivalent to the bipartite von Neumann entropy of the 
CQCP. 

The relationship between the CQCP entanglement 
spectrum, spin chain ground state probabilities, and CFT 
partition function is summarized by 
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Here T is the transfer matrix of the CFT, and the re- 
sult is independent of the boundary state \i) in the limit 
of an infinite strip. With a basis {cr} for the Shannon 
entropy, the projection operator \a) {a\ in the numerator 
of Eq. ([3]) constrains the corresponding field of the CFT 
along a cut at r = 0; the bulk fields remain free. The 
numerator and denominator are thus constrained and un- 
constrained partition functions, with free energies -F[{cr}] 



and F respectively. 

The Renyi entropies S'„ suggest a thermodynamic no- 
tation 

Z{n) = e-^(") = E^'- = E e-"(^»">l^^\ (4) 
{<y} {^} 

defined so that F{n) = (n — 1)5'„. The 'free energy' F{n) 
is found to depend on the length of the chain L as 



F(n, L) = h{n)L + 7(n) log(L/a) + /o(n) + 



(5) 



At a critical point ^{n) should be universal and deter- 
mined by the CFT as it is independent of the UV cutoff 
1/a. In the case of a cylinder geometry, 7(n) is zero for 
an arbitrary CFT as the trace anomaly from a smooth 
conformal defect line is zero. For a strip geometry, a uni- 
versal logarithmic term is expected as explained below. 

For integer n, Z(n) can be interpreted as the partition 
function of a replicated CFT with n copies of the field 
constrained to agree along r = 0, normalized by the free 
partition function. F(n) is then interpreted as the free 
energy of a defect line in the replicated CFT. If the loga- 
rithmic contributions can be calculated in the replicated 
CFT and analytically continued to n = 1, then the entire 
spectrum of Renyi entropies Sn as well as S is known, via 
S — dnF{n)\n=i. (Below, this 'replica trick' is found to 
fail for the Ising model due to non-analyticity of F{n) at 
n = 1 in the thermodynamic limit.) 

We start with numerical evidence for the existence of 
logarithmic terms in the case of the c = 1 compact boson, 
the CFT associated with the CQCP of the square-lattice 
quantum dimer model. For a free boson, Ref. |T3 pre- 
dicted a contribution —j log(L) S 5 in the strip geome- 
try; we indeed find a logarithm — ^ log(L) for the external 
boundary conditions (on the edge of the strip) assumed 
there, and for an alternate boundary condition we obtain 
+ \og{L) as derived below. By using the transfer matrix 
mapping, it is sufficient to calculate the Shannon entropy 
of a ID quantum model, the spin-1/2 XXZ chain 
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In the continuum limit the XXZ model is described by a 
compact boson. The boundary conditions are tuned by 
an external field h on the surface sites, and the Shannon 
entropy is calculated using the cr^ basis. 

The ground state was found by TEBD for L = 
{2,4, • • • ,30} at A = {-^-O, \, \, |}. We estimate the 
ground state error to be of order e = 1 — | (0 | 5*) p ^ 
0(10^®), which was verified for the exactly solved A — 
case. The partition function Z{n) was then computed 
by summing over all 2^ configurations. The logarithmic 
contribution 7(n) was extracted from the finite-size data 
as described in Fig. [2^. 7(n) is well described by 
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FIG. 2. a) The logarithmic part of F{n) at A = -1/2. Data 
is included both without (h = 0) and with {h > 0) an edge 
magnetic field. F{n, L) was computed for L = 2, 4, • ■ ■ 30. 
Windows of 6 data points centered ai £, L — {£ — 5, ■ ■ ■ , ^+5}, 
were fitted to the form /iL + 7log(L) + /o + /-ii~^. The data 
displayed includes 7 for £ — {9, 23}. The approach of 7(n) to 
the dashed lines suggests convergence to 7(n) = =p|(n — 1) 
as f — >■ 00, for h = and h > respectively. The inclusion 
of a term accelerates convergence, but the value of 7(n) 
remains insensitive to within a few percent to the choice of 
fitting form for terms that vanish as L — > 00, as discussed in 
the supplementary material, b) The logarithmic contribution 
to the von Neumann entropy S{L), as a function of the fitting 
window center £. As £ increases, the logarithm's coefficient 
converges to Tl/4, depending on the boundary conditions. 
Successive lines are for A = { — 1/2, 0, 1/2}. 



for h ^ and h > boundary conditions respectively. 
For the compact free boson, Eq. ([t]) is correct for all n, 
while for XXZ and dimer lattice realizations, the result is 
modified in the regime n > 2 due to additional boundary 
operators becoming relevant [T8| . 

The slope of 7(71) at n = 1 is of special interest as the 
von Neumann entropy is given by 5 = dnF{n)\n=i if the 
derivative exists. In the case of the XXZ model, there is 
a logarithmic contribution 9„7(n)|„=i = =F^, so the von 
Neumann entropy scales as 



S{L) ^s,-Lt\ log(i/a) 



(8) 

with increasing 



The rapid convergence of 5(L) to Eq. 
L is illustrated in Fig. [2]d. 

To derive Eq. ([7|, we perform a Jordan- Wigner trans- 
formation on the XXZ spin chain and bosonize the re- 
sulting fermionic model. The model is mapped onto the 
universality class of a free, compact boson </> for |A| < 1, 
with a compactification radius R that depends on A. For 
h = the external boundary condition (b.c.) is Dirich- 
let in the continuum limit, (j) — 0, while for h ^ the 
external b.c. is Neumann, dxcf) = [19]. The resulting 
change in the logarithmic contribution is thus attributed 
to 'boundary condition changing operators' (bcc opera- 
tors) in the underlying CFT. For integer n, the repli- 
cated theory contains n copies of the field, 0", subject 



to the constraint 0° = 4^ along the r = cut, modulo 
compactification. In the analogous problem for a non- 
compact boson, the action, external b.c.'s, and path in- 
tegral measure are all invariant under an 0{n) rotation 
of the replica index, so we can rotate the replicated fields 
to a new basis which includes the 'center of mass' field 
4>CM — ^X]a=i'^'^- The gluing condition now factor- 
izes, giving the free center of mass field plus (n — 1) de- 
coupled fields satisfying Dirichlet b.c. onthecut[T2l The 
free energy is 



F{n) = {n-1){FD~F) 



(9) 



where Fjj is the free energy of a single replica with Dirich- 
let b.c. on the cut and F is the unconstrained free energy. 
Using this rotation, the logarithmic contribution to F{n) 
follows from the free energies Fjj and of a single replica. 

Logarithmic contributions to the free energy of a CFT 
arise from a 'trace anomaly,' in which the trace of the 
stress tensor (T^) does not vanish. In the relevant 'cut 
strip' geometry, the anomaly arises from the four corners 
where the cut along t — terminates into the exter- 
nal boundary, as illustrated in Fig. [T^. The geometry of 
the corner contributes a term proportional to the central 
charge c [20]. However, there can also be a contribution 
due to the changing b.c: if the external b.c. is Neumann, 
then termination of the Dirichlet b.c. on the cut into the 
external Neumann b.c. introduces a bcc operator [21j. 
The scaling dimension of this bcc operator, hjsiD, will 
appear as an extra contribution to the trace anomaly. 

For an arbitrary CFT in the 'cut strip' geometry the 
predicted coefficient of the logarithmic term is 



LdLiF, - F) = 2 X 2 



(10) 



F is the free energy of the strip with no constraint at 
T — 0, while Fc is the free energy when constrained to 
some b.c. 'c' along the cut, and both obey external b.c. 
'e'. Here h^c is the scaling dimension of the bcc operator 
required to take the external b.c. e to the b.c. c on the 
cut. This result can be extended if e.g. the cut is not 
perpendicular to the boundary. In the compact boson 
case, c = 1, and the scaling dimension h^D is or 1/16 
for e — D,N respectively. Substituting into Eq. (10 1, 
LdL{F£) ~ F) — Tj for D/N respectively, i.e., the sign 
of the logarithmic term depends on the external b.c. The 
XXZ results are well described by the rotation trick and 
this bcc operator effect. It can be shown (supplementary 
material) that the logarithm is unaffected by compact- 
ification, as supported by the XXZ data. Briefly, field 
configurations can be decomposed into two sectors: fluc- 
tuations versus zero modes or vortices. Only the latter, 
topological sector is affected by compactification. The 
fluctuating sector can be rotated. The topological sec- 
tor has exact conformal invariance: the trace anomaly 
and logarithm arise only from the fluctuations. Hence 
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FIG. 3. The logarithmic part of F(n) for the TFI model. 
F{n, L) was computed for L — 2,4, - ■ ■ 40, and 7(n) extracted 
by the same procedure as for the XXZ model. The dashed 
lines are 7(n) = — |n and 7(n) = |n. In the inset, the data 
is collapsed by plotting 7(n)/n as a function of (n — 1)7'(1). 



the logarithmic terms are compactification-independent 
while the 0(1) terms, arising from both sectors, are not. 

The c — 1/2 Ising model exhibits strikingly different 
behavior from the compact boson. In the transfer matrix 
formalism, it is sufficient to consider the Shannon entropy 
of the critical transverse- field Ising model. 
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The Shannon entropy is calculated in the cr^ basis. 

The extracted logarithm 7(n), shown in Fig. |3j sug- 
gests a discontinuous n-dependence at n = 1. The dashed 
lines illustrate the likely convergence to 



7(n) ^ 



for n = 1 
for n > 1 



(12) 



Note that the rotation used in the boson case predicts 
that 7(n) c>c (n — 1), as there is always one 'free' center 
of mass field. In the present Ising case, "f{n) oc n, which 
is evidence that the rotation trick is not applicable and a 
new analysis is required. In the inset of Fig.[3j the data is 
collapsed by plotting ^{n)/n as a function of 7'(l)(n— 1). 
Recall 7(1) — identically due to the normalization of 
the entanglement spectrum. 

The large and small n limits can be understood as 
follows. At large n, Z{n) is dominated by the probabil- 
ity = I (tt • • • I 0) P and its partner p^. In this 'low 
temperature' regime, the defect line at r = is asymp- 
totically a 'fixed' defect line. Now all n fields experience 
this fixed defect condition, not n — 1 of them, and 



F{n) ^ n{F^ - F) - log(2) 



(13) 



Here F^ is the free energy of a single replica with spins 
constrained to cr = -1-1 along the defect. Again, there is 



an anomaly due to the 'cut strip' geometry. The exter- 
nal boundary conditions are free, so the relevant scaling 
dimension is hf^ = [21]. Using Eq. (10), 7(n) = |n 
as observed. Likewise, in the small n limit, the defect 
becomes disordered. This is obvious at the n = 1/2 
point, where the geometry is in fact a half-strip with 
a free boundary at r = 0. With hf f — 0, we arrive at 
7(n) = — |n as observed. While these results strictly ap- 
ply only in large and small n limits, the data appears to 
support an n-dependent phase transition: for n < 1, the 
defect flows to free boundary conditions for 2n decoupled 
half strips, while for n > 1, the defect is fixed. This in- 
terpretation is consistent with numerical results for the 
0(1) term on a cylinder ^l6i . The renormalization group 
equations for these b.c. flows can be analyzed perturba- 
tively [22], and support the hypothesis of a phase tran- 
sition. This perturbative approach could be generalized 
to other CFTs in order to study further possible entan- 
glement boundary conditions and replica transitions. 
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FIG. 4. The extracted logarithm 7'(n = 1) as a function of fitting window, for three fitting procedures: a) No subleading 
terms, b) a subleading term /_iL~^, and c) subleading terms f-iL^^ + f-2L~^. The data is taken from the XXZ model at 
A — J with Dirichlet boundary conditions. 

Supplementary material for "Logarithmic terms in entanglement entropies of 2D 
quantum critical points and Shannon entropies of spin chains" 

Fitting Procedure 

In order to extract the logarithmic part of F{n,L) = fi{n)L + 7(n)log(L) + /o + • ■ ■ from finite size samples, 
L = {2, 4, • • • }, we performed a windowed fitting procedure: each window of 6 data points, L = — 5, ^ — 3, • • • , ^ + 5}, 
was fit to the form fi{£,n)L + 7(71, £) log(L) + fo{n,£) + f^i{n,i)L~^ using linear least squares. As i increased, the 
coefficient ^{n,£) of the XXZ model was observed to converge for the range of n studied, 0.25 < n < 2, as shown in 
Fig. 2 of the main text. Here we address the sensitivity of the result to the details of the fitting procedure. 

In principle, it is not necessary to include further subleading terms such as L^^ if large enough system sizes are 
considered, but convergence can be accelerated by inclusion of terms arising from irrelevant operators. The choice 
of is motivated by the scaling dimension of the stress tensor, which we generically expect to be present. To 
demonstrate the improved convergence, we perform the windowed fit with a) no terms beyond /o, with a window size 
of 5, b) a subleading term f-iL~^, with a window size of 6 and c) subleading terms f-iL~^ + /_2i~^, with a window 
size of 7. We extract j'{n = 1,£) using all three procedures at A = ^, with Dirichlet boundary conditions, as shown 
in Fig. [ij For the largest available window (with an upper bound of L = 30), the three fitting procedures produce 
7'(n = l,£max) = {—0.233, —0.247, —0.249} for cases a), b), and c) respectively. In all three cases, the trajectory is 
consistent with convergence to 7'(n = 1) = —0.25 ±0.01. While including L~-^ and its descendants clearly accelerates 
the convergence, there are a multitude of other irrelevant operators, such as the vertex operators, with n - dependent 
scaling dimensions. To test the sensitivity of the fit to the presence of different power law type terms, we replace 
L^^ with a term L^°', with a 7 point window. We find that for a € (0.5,3), 7'(ri — 1,£ — 25) ranges between -0.248 
and -0.256, with trajectories again consistent with j'{n = 1) = 0.25 ± 0.01. To an accuracy of several percent, our 
conclusion is insensitive to the choice of subleading terms. 

A similar analysis holds for the Ising model, except near n = 1, where the extracted logarithm is not yet close to 
converging for the system sizes considered. In the approximate region 0.8 < n < 1.25, the extracted logarithm is 
sensitive to the choice of subleading terms or the introduction of terms such as log(L)^. This is to be expected in 
any region far from convergence. However, a monotonic approach to the extrapolated values 7(n) = {— |n, |n} with 
increasing £ is a generic feature, as well as a slowly diverging derivative 7'(1) at n = 1. 

Renyi Entropy of the Free Compact Boson CQCP 

The method applied by Ref. [T]to the non-compact boson does not strictly apply in the compact case. Several works 
have addressed the computation of the 0(1) term in the case of an infinite cylinder using the boundary conformal 
field theory of the compact boson. A recent preprint [S^ computes the full partition function, again in the infinite 
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cylinder or strip geometry, without using the repHca trick, with the advantage of applying immediately to non-integer 
n. Here we evaluate the replicated partition function for an arbitrary geometry. We show that the logarithmic part is 
not changed by compactification, and give a formula for the 0(1) contribution that can be extended to any geometry. 
Only for finite size subsystems is the effect of compactification really present, in the form of relative winding numbers 
which lead to Riemann Theta functions. In the limit of an infinite cylinder, we recover the formula for the 0(1) term 
found in previous work. 

We must evaluate the partition function of the replicated, constrained boson 

Z(n)=Z-"/" P[0'^]e^^J=i^[^°l (1) 

J con 

for a compact boson of radius R and action 

S[^] = ij,jd'x {V<pf . (2) 

We assume that the entanglement cut is a connected curve 7; the disconnected case (for a siibsystem of two disks, 
for example, or an annulus) can be solved in a similar manner. In the connected case, the constraint takes the form 
(/)°(a;) = (f)^{x) + 2'KRw"^ along the cut, where w"'^ G Z are a set of relative winding numbers. There is necessarily some 
ambiguity in the normalization of this constraint, leading to an extensive non-universal contribution F{n) ^ ^(n — 1), 
which, for our purposes, can be ignored. 

We assume an arbitrary mix of Dirichlet and Neumann external boundary conditions, suitably generalized to the 
compact case. The field configurations can be decomposed according to — (f)^^ + 0", where V^(/)q, = encodes the 
vortices and zero modes and is a fluctuation. Both obey the external boundary conditions. The fluctuations are 
not compactified. Compactification affects only the topological sector, with a indexing whatever winding numbers or 
zero modes are required. 

The action S[(j)] = S[4>a] + S[(j)] decouples into the topological and fluctuating part, so the replicated partition 
function becomes 



Z(n) 



E 



(3) 



The constraint imposes itself both in the topological and fluctuating sectors. For a fixed topological configuration a"", 
the fluctuations arc constrained along the cut according to 4>°'{x) — 4'^{x) = (j}^b{x) — 4>a.a{x) + 2i:Rw"''^ . The relative 
winding numbers w must also be summed over. 

For fixed a° and w"-^ , we have a linear constraint on the fluctuations of the form M ■ (/){x) — c{x) with less than full 
rank: the 'center of mass' mode J2a There are many possible choices of M and c; for example, we can 

choose a cyclic form 



M- 



/-II 
-1 1 
0-1 

V 1 



-1/ 



/ (j)ai (x) — 4>o,2 (x) + 2ttRw^ \ 



\ (j)an {x) - (pa^ (x) + 2-kRw'^ J 



(4) 



with relative winding numbers w G : (1, 1, • • • , 1) • w = 0. For any choice of M, M has a right null vector 
(1, 1, • • • , 1) corresponding to the 'center of mass mode', and a left null vector we call d. For the above choice of 
M, d = (1,1,-- - ,1). Note that while a particular choice of M may appear to break the permutation symmetry, the 
imposed constraint does not. 

As the fluctuations are not compactifled, we can choose some 0{n) rotation cf)' = Tcp such that ^ ^ is the free 'center 
of mass' mode. The measure and action 'D[^"']e~ «'je-I2S[4> "] are left unchanged, but the constraint 

becomes 



M ■ if){x) = c(x) ^ TMT-i • 0' (x) = T • c{x) 



(5) 



3 



By our choice of T, the center of mass ^ does not appear in the constraint condition, so wih contribute the 
partition function of unconstrained fluctuations which we denote by Zf. The remaining (n— 1) fields are constrained 
according to {x) — TM~^c[x) = c'{x), as M can be inverted in the subspace spanned by c. 

The precise form of c'{x) is irrelevant, except that it is given by a fixed linear combination of the original c{x). 
There is in addition a non-universal contribution ^/n to the partition function arising from careful consideration 
of the Jacobian of T acting on the original constraint: 



a" .w.con 



Again, we will ignore the non-universal extensive contribution. Finally, we use the shift invariance of the path 
integral to reduce the problem to Dirichlet constraints. Let C"''(x), defined on the entire region, satisfy the external 
boundary conditions as well as C'{x)\^ = c'{x) along the cut. Away from the cut we require V^C(a;) = 0, which 
uniquely specifies the C°. We now shift cj)' —i' (j)' + C . The redefined fluctuations obey Dirichlet conditions along the 
cut so the action again decouples, 

S[4>' + C] = S[4>'] + S[C'] (7) 



The partition function becomes 

aa,w:con 



The (n — 1) fluctuations decouple entirely from the topological sector, and are subject to Dirichlet constraints on the 
cut as wefl as the external boundary conditions. Let Zjj denote the partition function of the bosons obeying Dirichlet 
conditions on the cut, so that the fluctuations contribute -Z^^^. Furthermore, the single replica partition function 
also decouples into a fluctuating and topological component, Z — Zf ■ Ztop- The full partition function is 



Zin)^f^) ' ^ Y: 5:e-^--[^"l-^^;-«[^'] (9) 

The crucial observation is that S[C'] and ^[^q] are scale invariant and the sums require no UV regulation, so have 
exact conformal invariance. Any L dependent terms of Z„ necessarily come from the fluctuations, ( ^ ) j which 

are identical to those of a non-compact boson. Thus, as far as the log(L) term is concerned, the rotation trick gives 
the correct result. 

To compute the 0(1) contributions, we must calculate ^[C]. Define C{x)\~^ = c{x) in the same manner as the C'{x), 
which will be related according to C'{x) = TM^^C{x) (again, M can be inverted in this subspace). All dependence 
on T drops out, and we find 



n « 

Y^S[C'] = f- / VC^{x){MM^y^VC{x) (10) 

a=2 ■' 



With knowledge of the classical solutions C{x) for a given sector a, w, the 0(1) term can then be computed. 



0(1) Terms for the Cylinder Geometry 

We now reproduce the 0(1) contribution \og{^/2gR) — ^ [3l IH [6] on an infinite cylinder. We consider a cylinder of 
circumference L and length 2/3, with the entanglement cut at r = 0. We assume Dirichlet conditions dX t — ±/3, so 
that there are no topological indices a to sum over. 

A simple representation of the constraint is 
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/ 0\ 

-1 1 ••• 

-10 1 

: ■■■ 

V -1 I J 



( 2ttRw^ \ 



(11) 



Obviously = 0, corresponding to the left null vector d = (1, 0, • • • , 0) of M. It is simple to find the desired 
classical solutions, by choosing C{t,x) — c(l — ^-^), with c — 2ttRw. We find 



-^^(2nRf I w'^iMM'^y^w = ^(2TTRf^w^(MM^)-^w 
An p J 2tt p 



(12) 



Inserting into our general result, we have 



Z(n) 



E 



^-^^{2TiB.fw'^iMM'^)-'^w 



(13) 



w£7a^ -.d-w—O 



The sum is a Riemann Theta function, which we will evaluate in the large (3 limit. The lattice sum goes over to an 
integral 



w^Zi'^-.d.'w—O 



..n-l 



(14) 



where dw" ^ is the Euclidean measure in the subspace d.w — 0. Standard gaussian integral identities give 



Z(n) 



2f 



(«-l)/2 



|MM^|' 



\L2gR^^ 

where IMM'^I' is the determinant of MM'^ with the zero mode omitted. We find IMM'^I' = n. So 



(15) 



Zin) 



2f 



- o\ ("-l)/2 



(16) 



The asymptotics of Zjj / y (which are given by rj functions) show that the fluctuating part does not contribute at 
0(1). Furthermore, dependence on the modular parameter ^ contributes an overall divergent constant, so can be 
ignored. We extract the universal part 



Zin) - (v/2gi?)-("-i), 
with resultant entanglement entropy S = \og{y/2gR) — ^ as desired. 



(17) 
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